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Abstract 

Lattice measurements of spatial correlation functions of the operators FF and FF in thermal 
SU(3) gauge theory have revealed a clear difference between the two channels at "intermedi- 
ate" distances, x ~ 1/nT. This is at odds with the AdS/CFT limit which predicts the results 
to coincide. On the other hand, an OPE analysis at short distances (x <C 1/irT) as well as 
effective theory methods at long distances (x ^> 1/vrT) suggest differences. Here we study 
the situation at intermediate distances by determining the time-averaged spatial correlators 
through a 2-loop computation. We do find unequal results, however the numerical disparity 
is small. Apart from theoretical issues, a future comparison of our results with time-averaged 
lattice measurements might also be of phenomenological interest in that understanding the 
convergence of the weak-coupling series at intermediate distances may bear on studies of the 
thermal broadening of heavy quarkonium resonances. 
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1. Introduction 



Spatial correlation functions of local gauge-invariant operators (such as components of the 
energy-momentum tensor or various charge densities) offer a set of observables that can yield 
a detailed characterization of the physical properties of a thermal system. Consequently, 
in lattice QCD, correlators related to gluonic (flavour-singlet) as well as hadronic (flavour 
non-singlet) operators have been studied since dawn [TJ [2] and continue to be important in 
the present era of unquenched simulations (see, e.g., refs. 011]). 

At high temperatures (T 3> 200 MeV), spatial correlation functions can be sensitive to dif- 
ferent types of physics at different distance scales. At short distances (x <C 1/irT), we expect 
them to display vacuum- like behaviour; at intermediate distances (x ~ 1/ttT), we expect 
a complicated functional form with thermal modifications of order unity; at long distances 
(x ^> 1/irT), the correlators should ultimately reduce to a single exponential, characterized 
by collective phenomena such as colour-electric or colour-magnetic screening [5]. In addition, 
the zero-temperature confinement scale may also make an appearance, if a measurement at a 
low temperature is subtracted from the high-temperature measurement in an effort to remove 
the vacuum-like behaviour at short distances. 

For future reference, we note that the distance scale that one is interested in may influence 
some technical details of the lattice measurement carried out. For instance, for short-distance 
physics a signal may be obtained from point-point correlators, whereas the long-distance 
asymptotics can perhaps best be extracted by using operators averaged over the Euclidean 
time direction as well as over the transverse spatial plane. 

In the present paper, following a recent high-precision lattice measurement [6], we consider 
spatial correlators within pure SU(3) gauge theory at a temperature above that of the decon- 
finement transition. The operators correlated are Tr [F^F^] (corresponding physically to 
the trace of the energy-momentum tensor, modulo a coefficient) and Tr [F^F^] (correspond- 
ing physically to topological charge density, modulo a coefficient). We assume the operators 
to be averaged over the Euclidean time coordinate but to remain point-like in space coordi- 
nates. On this point we actually differ from ref. [6] in which the operators were point-like 
also in time; apart from technical simplifications, our choice is motivated by the fact that 
the correlators then become more analogous to Polyakov-loop ones which are relevant for 
the phenomenology of heavy quarkonium physics (the ramifications in this direction will be 
discussed in the conclusions). Time-averaged correlators are also the ones containing the 
thermodynamic information relevant for spectral sum rules (see, e.g., refs. 0|8]). 

One of the main findings of ref. [6] was that, once the (genuine) vacuum part was subtracted 
from the correlators, the two channels showed strikingly different behaviours at intermediate 
distances (cf. fig. 5 of ref. [6J): for Tr [F^ U F^ U ] the thermal modification was negative in 
this range (reducing the correlator) whereas for Tr [F^F^] it was positive (enhancing the 
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correlator). That such a difference exists is in contrast with the strongly-coupled AdS/CFT 
limit of J\f = 4 Super- Yang-Mills theory which predicts the two correlators to coincide [9j[6]. 

Within pure SU(3) theory, the emergence of a qualitative difference between the two chan- 
nels has previously been understood both in the short-distance (x <C 1/nT) and in the long- 
distance (x S> 1/nT) limits. In the short-distance limit, the correlators can be expanded in 
an Operator Product Expansion (OPE). The leading terms are vacuum-like contributions, 
which already differ as far as non-logarithmic terms are concerned [10] . Furthermore, OPE 
contributions proportional to the expectation value of the trace anomaly, (e — 3p)(T), come 
with opposite signs |11[ I12j. underlining that any breaking of conformal invariance is bound 
to lift the degeneracy. 

On the long-distance side (x 3> 1/nT), the very fact that the two operators have dif- 
ferent discrete quantum numbers has been exploited for defining colour-electric and colour- 
magnetic screening masses in a non-Abelian theory [5]. The colour-electric channel couples to 
Tr [FnvFftv] and the colour-magnetic to Tr [F^F^]. Practical measurements in pure Yang- 
Mills theory and QCD (making use of an infrared effective field theory) can be found in 
ref. [13] and indeed show substantial variation in the screening masses; the qualitative differ- 
ences in the corresponding weak-coupling expressions have been elaborated upon in ref. |14j . 

The purpose of the present paper is to "interpolate" between the known short-distance 
and long-distance limits, by considering intermediate distances x ~ 1/nT. This situation is 
"simple" in the sense that straightforward perturbation theory is formally applicable, without 
the need of re-organization in terms of OPE or of effective field theories. The challenge is 
that the correlators have a relatively complicated functional form. Nevertheless, we have 
managed to obtain expressions for both correlators, and foresee an opportunity for a practical 
comparison of the weak-coupling expansion against lattice simulations of the type in ref. [6]. 

The plan of this paper is the following. The observables considered are specified in sec. El 
and the main steps of the computation are described in some detail in secO Section [4] collects 
together the basic results applicable for short and intermediate distances, whereas in sec. [5] 
we elaborate on how the perturbative expressions need to be resummed in the long-distance 
limit. Numerical illustrations comprise sec. El and our conclusions are collected in sec. [7J 

2. Setup 

Our basic notation follows ref. |12j . so we discuss the setup only briefly. We consider pure 
SU(iV c ) Yang- Mills theory (with N c = 3 in numerical estimates), dimensionally regularized by 
analytically continuing to D = 4 — 2e space-time dimensions (Se = Jq^-t J d 3_2<E x j F^F^). 
The operators considered are 

= C 9 glF^F" , X = C X ^^vpa9i F av F pcr J C 2 - 1 ) 
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where g 2 is the bare gauge coupling squared; these operators require no renormalization at 
the order of our computation. We do not need to specify the coefficients c$, c x , but note that 
in ref. [6] the values eg ~ — ^ j c x = were chosen, where bo = ^n^y2 ■ The bare gauge 
coupling squared can be expanded in terms of the renormalized one, g 2 , as 

llg 2 N c 1 



1 



3(4vr) 



+ ... 



(2.2) 



where A denotes a scale parameter. The MS scheme renormalization scale is denoted by 
A 2 = 47rA 2 e _7E . In order to avoid unnecessary clutter all appearances of A 2e , which cancel 
when the result is expressed in terms of the renormalized gauge coupling, will be suppressed. 

We wish to evaluate the behaviour of the correlators of and x as functions of the spatial 
distance x = |x|, integrated over the Euclidean time r. To this end, we define the quantities 
{X = (r,x); P = (p n ,p); Gg(X) = (d(X)0(0)) c ; G X (X) = (x(X) X (0)» 



G (x) = / drGe(X) , 
J o 

Using the corresponding Fourier transforms 



G x (x) = 



drG x (X). 



Gg(P) 



-iP-X 



G e (X), G X (P) 



x 



iPX 



G X {X) , 



(2.3) 



(2.4) 



x 



we see that for our purposes it suffices to consider the static (p n = 0) limit of the correlators 
in momentum space: 



Gg(x) 



e i v*G e (pn = 0,p) 



G x {x) 



Jp 



e i P x G x (p n = 0,p) . 



(2.5) 



We work up to next-to-leading order (NLO) in perturbation theory. A natural starting 
point for the computation are then the expressions in eqs. (3.1) and (3.2) of ref. |12j . ob- 
tained after "scalarizing" the 1-loop and 2-loop Feynman graphs (fig. [T|) contributing to the 
momentum-space correlators (d,A = X 2 — 1): 



Go(P) 



(D-2) 



4dA<4sf: 
+ g 2 B N c {2(D-2) 



22 - ID 



- Ja + ^ 



-{D-l)Z & + {D-A)Z h 
(D - 4) 2 



G X (P) 



-16d A c 2 gi(D - 3) 
+ glN c {2(D-2 



2 

(D-2) 



-2g + (£>-2) 



+ {D-2) 2 l c -l d 
3X e + 3X h + 21i - 2j 



(2.6) 



-Ja + ^ <Jb 



-2a + (D- 4)2 b 



+ (D - 2) 2 



2L> 2 - 17L> + 42 



2f - 2(D - 4)2 g + (-D - 2) 



2 C — 2d 

32 e + 32 h + 22; 



(2.7) 
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(ii) (iii) (iv) 



(!) 



O oo <X> 



(v) (vi) (vii) 



Figure 1: The LO and NLO Feynman graphs contributing to the correlators of eqs. (|2.6p . (12.71) . 

The definitions of the master sum-integrals J & . . . Zj can be found in eqs. (|3.1f ) — (j3.13f) below. 

It will become clear later on that at long distances, the NLO corrections can overtake the 
leading order (LO) terms, indicating a breakdown of the perturbative series. In this situation 
perturbation theory needs to be resummed through effective field theory techniques. We 
return to when this happens in sec. [5J for the moment we simply compute the correlators as 
they stand in eqs. (j2.5[) — (j2.T|) . The results will be denoted by Gg , with "f" standing for 
"full theory". 

3. Naive full theory computation 

In this section, we describe in some detail the evaluation of the various master sum-integrals 
appearing in eqs. (|2.6p . (|2.7| . once they are Fourier transformed according to eq. (|2.5p . 

3.1. Sum-integrals 

The master sum-integrals are defined using dimensional regularization in D = 4 — 2e dimen- 
sions and, as in ref. [12], read $ Q = TY^ qn J q ) 




(3-1) 



(3.2) 



(3.3) 



(3.4) 



4 



Xc s 1^ (3 - 5) 

Id 55 t Q , R q 2 r\r - py ' (3 - 6) 

Xc = %. R Q 2 R 2 {Q-R) 2 ' (3 ' 7) 

Xf = ^ Q*(Q-R)*(R-P)* ' (3 ' 8) 





p)2 






Q 2 R 2 (Q - R) 2 (R - 


P) 2 


(Q - Pf 




Q 2 R 2 (Q-R) 2 (R- 


P) 2 


4(Q-P) 2 





Xh - %, R Q 2 R 2 (Q-R) 2 (R-P) 2 > (3 - 10) 
3 - i. ~,^1 Q ~ P L ^ , (3-11) 



Zi ' ~ %, R Q 2 R 2 {Q-R) 2 {R-P) 2 ' (3-12) 

Xj = t Q , R Q 2 R 2 {Q - R) 2 (Q - P) 2 (R- P) 2 ' (3 ' 13) 

As specified in eq. (|2.5p we are interested in taking their three-dimensional (3d) Fourier 
transforms, i.e. the (dimensionally regularized) functions 

l m (x) = [ e i v*l m (p n = 0,p) . (3.14) 

We adopt a notation in the following whereby p stands for a four-vector with vanishing 
frequency, i.e. p = (0, p); in the above definitions, we can thus replace P 2 — > p 2 , (Q — P) 2 — > 
(p - Q) 2 , and (R - P) 2 -> (p - P) 2 . 

In some cases, the contribution of the zero modes (q n = r n = 0) to the master sum-integrals 
is delicate, and a proper handling requires the use of a consistent regularization procedure. 
We therefore separate the zero-mode contribution from the rest, and evaluate it in sec. 13.21 
the results are denoted by 2^. In subsequent sections, we assume that the zero modes have 
been subtracted from the sum- integrals; the subtracted sum-intergals are denoted by 2^. 



3.2. Zero modes 

Setting q n = r n = 0, we see that many of the master sum-integrals vanish as scaleless integrals 
in dimensional regularization, 

Jl = Tl = l° h = l° c = Z d ° = l° e =0. (3.15) 
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For the rest, we first note that by dimensional analysis, 



J b ° ~ T(p 2 )-\ (3.16) 

l° m ~ T 2 {p 2 ) x ~ 2e , me {f,g,h,i,i',j}. (3.17) 

Recalling on the other hand that 

l ^ p 3~2e = V_ + Q{e) I e ip x p2 -4, = __6e + 0(e2) _ ^ 
Jp TT X° J p 7TX° 

we observe that it suffices to evaluate the master integrals 2^ to order 1/e. An elementary 
calculation (making use of the "zero-temperature parts" of the various master sum-integrals 
given in appendix A of ref. [12], after setting D — > 3 — 2e in them) now produces 

J° = 9 -^+0(e), (3.19) 
3tt 2 T 7 

2f° = — ^r-+°( e )> ( 3 - 2 °) 

X g ° = 0(e), (3.21) 
3tt 2 T 7 

X h ° = _^+0( e ), (3.22) 
3tt 2 T 7 

2° = 2ZL_L +0(e ), (3.23) 

2° = 0(e), (3.24) 

6tT 2 T 7 

2° = + 0(e), (3.25) 
where we have defined the dimensionless variable 

x = 2vrTx , (3.26) 
in terms of which we will present our results. 

3.3. J a and J h 

The 1-loop sum-integrals J & and Sh are only needed in the combination 

,2 
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where we have discarded a contact term. The remaining term can be evaluated in a straight- 
forward manner by writing p 4 as a spatial derivative and using results derived in ref. |15j . 
We will later on also need the 0(e) part but, for the time being, we merely write 

- - 4vr 4 T 7 d 4 /cothx\ , , . 

Jb _ 2J , = __(_j +0(£) . (3.28) 

Note that this expression contains both the zero-mode and the non-zero mode contributions; 
the zero-mode contribution of eq. (|3.19p is recovered as the only power-law term at long 
distances, if we expand cothx = 1 + 0(e~ 2x ). 



3.4. X a , 2T C , X e 

In momentum space, the P-dependence of the sum-integrals Z m , m £ {a, c, e}, factorizes from 
the Q and R dependences (in fact, I e vanishes identically). One is then left with taking a 
3d Fourier transform of a non- negative integer power of the momentum squared, which leads 
to a result proportional to the delta function or its derivatives. The results could be directly 
extracted from ref. p2] , but in the present paper we are only interested in finite values x ~ 1 , 
so we set these integrals to zero, 

Z m ^0, m€{a,c,e}. (3.29) 



Contact terms also appear as parts of the other master sum-integrals, and will be omitted 
there as well; again, if needed, they can be extracted from ref. |12j . 



3.5. Xb, Zd, Xg 

Next, consider the sum-integrals Z m , m S {b,d, g}. All of these consist of 1-loop integrals 
of the type that have been thoroughly investigated in the literature, see e.g. refs. [1"5]-[17|. 
Neglecting again contact terms, and noting that in eqs. (12.6p . (12.7P the sum-integrals Zb and 
Z g are multiplied by terms of 0(e), we obtain to the required order in e: 

Z b = 0(e°), (3.30) 
vr 2 T 7 1 d 3 



Xd = "ir^( C0th ") + 0(e) ' (3 - 31) 



It is worth noting that the contribution of the zero mode q n = r n = vanishes or is of 0(e) 
in all of these cases, so it does not need to be separately subtracted. 



2 T 7 1 d 4 /cothx \ 1 
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3.6. T f 



As a warm-up for the more challenging integrals, let us consider X[ in detail. Defining 

m)s imhw <3 - 33) 

and using a representation from ref. |T5] (cf. also eq. (|3.4ip below, with the flag a set to 
zero and the zero- mode part TS rn /8r which has already been accounted for in eq. (|3.20p 
omitted), we readily obtain for the sum-integral in question (the prime reminds of the zero- 
mode omission): 



T' 



-V 2 



-V 2 



n(g) 

r(p-R) 2 



3 i(p-r)-x 



p Tr(p-R)' 



(4vr) 2 \R 2 



£l 1^ + 2 



+ 



T 



(4vr) 5 



d 3 y 



e «r-(x+y) 

y 2 



coth y - - - 5 rr | < 
V 



\rn\y 



(4^) ; 



A-i + 2 £1 
m \ -> e 

+ T E- 

37T 2 r 7 id 2 / i 

4 x dx 2 V x 2 sinh 2 x 



x 

-2\r n \x 



{R 2 r 



( coth x — — — 5 r 



r 



x 



(3.34) 



In the last step we expanded in e and dropped contact terms; we also inserted the integrals 



-|jn|x 



r 7-2 + r 2 



47TX 



e ir x ln(r 2 + r 2 ) 



1 



27TX 3 



(1 + \r n \x)e 



\r n \x 



and carried out subsequently the elementary sums over r„ 



y^e 2|r ' l|x = cothx , ^|r n |xe 2|r " |:r 



x 



2 sinh x 



(3.35) 



(3.36) 



3.7. 

Next, we inspect the sum-integral 1^, the evaluation of which is in many ways similar to that 
of If. Separating again the zero temperature part of Il(i?) from the rest, we obtain 



T' 



A 



2c 



(4vr) 5 



+ 



(Air) 



1 

- + 2 



rV 4 



p 

jpx 



rR 2 (p-RY 



dV 



colli// - <L ) e~l''" ;/ 



fO ,fT 



(3.37) 
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where we leave the first part intact for now (we return to it in sec. 13. lip and only consider 
the second term, TF, For this function, we obtain 



fT 



T 2 V 4 



d 3 y 



V |x - y| 



coth y — — 



IT T d 4 



2x dx 4 Jo xy [ V y 



y 

y< + ln 



e 2y> _ i 
e 2(x+y) _ i 



y< 



(3.38) 



where y < = min(y, x) and y> = max(y, x). After taking the derivatives with respect to x, we 
are left with an integral that appears to be difficult to perform analytically, but is numerically 
quite benign. Thus, we have opted to evaluate it numerically (cf. sec. H]). 

3.8. X[ 

As in ref. [12] , we reduce the evaluation of X[ to the other integrals, in particular I[, , employing 
the relation 

i! = i' a + i' e -%+i!,. (3.39) 

3.9. 1[, 

The evaluation of T[, involves a new type of 1-loop "self-energy" structure, as instead of H(R) 
of eq. (j3.33p it contains the function 



iWi?) EE 



qQHQ-RY* 



A(R) 5ij + B(R) n r 



(3.40) 



where in the latter step we have used three-dimensional rotational symmetry to write the 
tensor in terms of two scalar quantities. As usual A and B can be determined by projecting 
the defining equation with dij and r^r-,-; moreover, by completing squares, the results can be 
reduced to the known tadpole f, Q l/Q 2 = T 2 /12 as well as to the self-energies (cf. refs. [1~5|[1~T]) 



-qQ^q-ry 



(4vr) 2 \R?) \e ) 6R 2 8r 



2 _ r 2 



, d 3 v ■ 

1 /A 2 ^ e 



1 ay 

( coth y - d Tr 

y 3 



y 2 V 



■\i~n\y 



(3.41) 



Q 2 (Q - Rf 



T \ oTV 2 
(4vr) 2 V R' 2 J D - I V e + J + 6R 2 



(47T) 



d 3 y 



: d 2 

„2 V 



, _ 1 ay 
coth y — — 

y 3 



■\rn\y 



. (3.42) 



The coefficient a serves as a "flag", to be chosen freely (below a = 1 for r n ^ 0). Note that 
in eq. (|3.42p the zero- mode subtraction has actually no influence; 9y[l] =0. 
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Defining now, for brevity, 
fn(y) = 
fn(y) = 



^cothy - ^ - | - S rn \ C 



\rn\y 



cothy 



l _ y 

y 3 



1 



v 



5^ e"l' 



\r n \y 



(3.43) 
(3.44) 



where in the latter case we have, for future convenience, removed the a-flagged contribution 
in the case of the zero-mode, the functions A and B become 



A(R) 



R 2 



1 /A 2 \7l 



4(D-1) (4vr) 2 \R 2 J \e 



+ 2 



Tr 
64 



I n 



R 2 



B(R) 



g r 2 ( 47] -)2 



d 3 y^a 2 + r 2 )/ n (y) 



(3.45) 



1 r£vri + 2i + 



4(D-1) (An) 2 V«V V e 



1 _Srn.Ml + !E_ 



12 16 / i? 2 \ 64r 



+ 



1 T 

8T 1 (An) 2 



(3i? 2 - 2r 2 )c> 2 + (i? 2 + 2r 2 )r 2 f n (y) . (3.46) 



Defining A' = A + {^5 Tn } and B' = B — {^S Tn }, i.e. dropping the terms whose effects were 
already accounted for in eq. (I3.24p . the sum-integral under study becomes 



1 / e ip x 
'p 



1 



-4V 2 Tj^ 



rR 2 (p-R) 2 

e ~\r n \x 



Airx 



p 2 A'(R) + (p • r) 2 B'(R) 
K A>(R) 



R 2 



ipx 



4t^ -v 2 



,— |r n |a; 
Attx 



+ 



p p z + r, 



d 2 



■(r 2 + p-r) 2 B'(R) 



2 

n J r 



i? 2 



: ^CR) 
ii 2 / ' dxidxj 



2 



4irx dx-idxj 



B\R) 
R 2 



(3.47) 



which defines the two functions we now set out to compute. 

Let us start with and divide it into two parts according to the two terms in eq. (|3.45p . 

A — n-AO | n-AT 



where the first "zero-temperature" piece reads 



(3.48) 



rAO 



1 



A 2e /l 



4(D - 1) (4tt) 3 Ve 



+ 2] T^V 2 



6(4tt) 



V 2 



l^(l + jr n |x) e" 2 ^^ 



7T 2 T 7 1 d 3 /coth: 



48 a; dx 3 \ x 2 



(3.49) 
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In the second step we expanded in e, dropped contact terms, and made use of eq. (|3.35p ; in 
the last step the sums were carried out like in eq. (|3.36[) . The finite temperature part can, 
on the other hand, be written as 



r-AT 



T 



8(4tt) 3 

T 2 
8(4tt) 3 

T 2 
8(4^)3 

T 2 



v 2 E- 

r„ 

2r e 



X 

-\r n \x 



4(4^)3 V ^ 



7T 2 T 7 d 2 

8x dx 2 



x 

y 3 



fn(x) 



X z 
fn(x) 



+ /d 3 y 



1 



4ny 2 \x + y 



dlUy) 



dy — d y fn(y) 
y> y 



coth 



, f n (y) 
dy — o— 

yd 

(coth 



V 2 



- 1 V 

y 3 



l 



(3.50) 



where the last integral again seems to require numerical evaluation. 

Moving then on to iff, we once more split the integral into two parts, now according to 
the T = and terms of B'(Q). This produces 



rBO 



i*° + i{ 

D 



BT 



K 2e 



4(D - 1) (4tt) 2 \e 



+ 2 



jpx 



(P • r) 2 



(3.51) 



which we analyse together with the other divergent terms in sec. [37TT1 With Jj? T , on the 
other hand, we encounter a similar albeit more complicated expression as with IjF. Its 
evaluation nevertheless proceeds along the same lines: One begins by performing the spatial 
r integral, then takes care of the angular part of the y integration, partially integrates, and 
finally performs the sum over r n (which we leave intact for the moment). The result of the 
procedure reads (a prime indicating the omission of the zero-mode) 



r-BT 



Tw d 2 e-\ r "\ x d 2 
4^2^ 



dxjdxj 



x 

"OO 



rp2 e -\r„\x 



dxidxn I 12 8ir\r„ 



8(4vr) 



2 f° d tffl_U 1,22^ 
n\X Jo r 



-\r n \(x+y) 



fn(y) 



+ 



8(4vr) 2 \rl\x J y 3 



\rn\y + r 2 nV 2 ) e-^-^Uy) 



+ 



8(4tt) 2 

t 2 d 2 



°°dy 

y 3 
l 



y 
o 2 



X 



r* 2 I /v* 3 I ry* 

nr> * On « 



[1+ |r n |y + r 2 y 2 ) e" 



■knKy-a;) 



8(4-7r) 3 dxidxj x dxjdxj \ 3x ^ / 



d2/ 'j-^) /o(y) 



fn{y) 

. (3.52) 
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The derivatives can be taken by noting that the independence of the above expressions on 
angular variables implies 



didj g(x)didj f (x) 

Then the first row can be simplified into 

Q2 ( e -\r n \x Q2 



dl (* V") - (<?/') • 



.2 z 



(3.53) 



-y' 



dxidxj 

TT 2 T 7 



x dxjdxj 



rp2 e -\r n \x 

12 8vr|r n | 



n(2x) n'(2x) n"(2x) 
2x 4 x s x 2 



n^{2x) 



(3.54) 



where n stands for the bosonic distribution function n(x) = l/(e x — 1), whereas the last row 
(the contribution from r n = 0, q n ^ 0) becomes 
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8(47r) 2 dxidxj \ 4ttx dxidxj 



tt 2 T 7 
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+ 



(3.55) 



The other parts contain integrals that we collect together in sec. H]for numerical evaluation. 



3.10. Xj 

From the numerics point of view, the sum-integral 1'^ appears to be the most non-trivial 
one, as there are not many analytic tricks that we can apply to it. The most important 
observation is that 2j' is both UV and IR finite, implying that we can set e = and consider 
its evaluation in three-dimensional coordinate space. Fourier transforming all the propagators 
and subsequently rescaling the radial coordinates by x, we obtain 

V*' R -~\Qn\(y+z)-\r n \(\x-y\ + \x-z\)-\q n -r n \\y-z\ 



T 2 V 6 



d^y d A z 



yz 



(47r) 5 j yz\x — y||x — z||y 

~ [ dydz [ diidi 2 / d r 
16x J J_x J |e - y||e - z||y - z\ 

v d6 /™2V^' -s[|m|(jf+5)+|n|(|e-y| + |e-z|)+|m-n||y-z|] 



(3.56) 



where e is a unit vector and the prime on the sum means that the simultaneous zero mode 
m = n = has been left out. This sum can now be explicitly carried out, giving 

ttT 7 

I dydz I 



8.T 



/"OO /*1 /*27T 

/ dydz I dti dt 2 \ d 
'o J-l Jo 



y\\ e - z \\y 



(3.57) 



d 6 



x i—rS x 
dx 6 I 



n(xa) + n(xb) + n(xc) + n(xa)n(xb) + n(xb)n(xc) + n(xc)n(xa) 
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where we have defined the combinations a = y + z + |e — y| + |e — z|, 6 = y + z + |y — z| 
and c = |e — y| + | e — z | + |y — z|. The sixth derivative of the expression on the last line 
contains a large number of terms, but its evaluation is straightforward to automatize using 
the relation n'(x) = — n(x)[l + n(x)]. The remaining 5d integral can be performed using 
numerical methods. 



3.11. Divergent terms 

Above, we have left two integrals unevaluated, namely the UV-divergent term ZjJ of eq. (J3.37P 
and of eq. (j3.51|) . These contributions are expected to partially cancel against the 
renormalization of the gauge coupling in connection with the 1-loop terms in eqs. (|2.6p and 
(12. 7h . Here, we evaluate the sum of these three functions, obtaining a UV-finite result. 

Inspecting the forms of Gg(x) and G x (x) in eqs. (j2.6j) and (j2.7|) . we observe that the only 
parts containing divergences are 



r<div ridiv 



Jb + 5B^c(3^h + 2x4X i ?°) 



4d A c 2 e -\U A c\{B - 3) 
Expanding g\ according to eq. (|2.2[) . we are left to consider the combination 

i r e i P .^ i r up 4 



(3.58) 



(4^) 2 7 p %R 2 {p-Rf\ 36 

+ (H(^)> 4 +^r^ 

\2 ^4 



, ,p4 1 f 4(p-^- P 4 

3(4vr)2 J p %R?{p-Rf\ 



2 , ^ 

Writing the scalar products (p ■ r) 2 in terms of — R) 2 , we can massage the result into a 
relatively simple form; in particular, apart from contact terms, the 1/e divergence disappears 
in the course of the r-integration. Upon denoting 

D„(x,r„) s /(l^e'-^, (3-60) 

5„ (l , r „) s /|L,,*Hp!> , (3. 61) 
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the remaining finite result becomes 



-Z-div — 



2T 



3(4vr) 
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The various functions that appear here have the forms 
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which at the end leads to the result 
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4. Basic results 



4.1. General expressions 



We are now ready to write down our final results for the full theory contributions to the 
correlators Gg(x) and G^(x), from eqs. (I2.6p . (I2.7P and (I3.39p . Omitting contact terms and 
collecting terms up to O(e ), we obtain 



01 



Ad A c 2 
r<F 



g A J h + 2g 6 N c 
g 4 J h + 2g 6 N c 



Idiv — 2Zd — 3Zf + 3Zh + 2I;> — Xj + 2eZg 



(4.1) 
(4.2) 



It is immediately clear that the two channels only differ through a single term, the last one 
in eq. (|4.2p . Inserting the result from eq. (|3.32p . the difference has a simple form; it is shown 
explicitly in eqs. (|4.6p . (|4.7p below and is displayed numerically in fig. [2j 
Adding up the various parts from sec. [3] and defining the functions 



n(x) 
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i = E e ~ na 



n=l 
oo 



n=l 



rr 



i(x) = ln(l — e 

oo 

Li3(e-) = ^ 



E 

n=l 



?l 



n=l 



71 



(4.3) 
(4.4) 



we end up with the expressions 
G F e (x) 



Ad A c 2 e 



5 4 4vt 4 T>lo(x) 
2- r r i 11 



+ g N c n z T 



4>lo(x) In 



e 7E Ax 
2vrT 



+ ^ c 27r 2 T 7 </- LO (x) . 



4d A ci 



+ (p^ho(x) + 01 (s) + 02 (*) ? , (4.5) 



(4.6) 



The function </>LO; which determines the leading-order behaviour, the scale dependence, as 
well as the difference between the two channels, is given by (cf. eq. (|3.28p ) 



0LO(X) = 3g + - —± 



24 2 d 4 (n{2x) 
a; 6 x dx 4 
( 120 



'(- 



24 



+ e 



-2x 



X < 1 

X » 1 



(4.7) 
(4.8) 
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The function </>nlo(x), which determines the long-distance asymptotics of the next-to- leading 
order correction, reads 



0NLo(x) 



4 24 
+ -=■ + 



i(2x) - ln(2x) 
1 



3x 4 
( 64 40 

4 | 24 881n(2x; | g /- 2 , 



x < 1 



, x > 1 



(4.9) 
(4.10) 



3x 4 x 5 x 6 

The remaining functions 0i and </> 2 are subdominant at long distances (4>i by exponential 
terms, 2 by powerlike terms), and take the respective forms 
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yz 



dx 6 \' 



|e — y||e — z||y — z| 

n{xa) + 2n(xb) + n(xc)n(xb) + 2n(xa)n(xb) 



(4.12) 



y + z + |e - y| + |e - z| , 6 = y + 2 + |y - z| 



a + 6 - 2(y + z) 



The integrals still remaining in the definitions of (p\ and 2 are carried out numerically. The 
analytically known functions 0lo> 0nlo are illustrated in fig. El the numerically determined 
functions (f>i , c/> 2 in fig. [3l 

Let us work out the asymptotic behaviours of 4>i and </> 2 . For 0i, we find 



824 58 
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i>f(x) 



0(e 



-2x\ 



3x 4 45x 3 
x > 1 . 



+ 01 r=s ] , x < 1 , 



(4.13) 
(4.14) 



16 




Figure 2: The basic functions denned in eqs. (|4.7p . (|4.9p . together with their asymptotic limits as 
given in eqs. (|4.8[) . (|4.10[) . The functions UV / IR indicate the nature of the limit, and the argument 
shows the order up to which terms have been included. 



For <p2, in turn, the short-distance behaviour is 

The infrared behaviour is, in contrast, rather non-trivial: while all terms in (fi\ are exponen- 
tially small at long distances, there are some parts in <p2 which are not. Indeed, despite a 
large x, there is a slice in the integration range of eq. (|4.12|) where b is small (because of 
the appearance of the unit vectors, this is not possible in the terms involving a or c). The 
intuitive meaning of this term is perhaps clearest if it is rewritten in momentum space: 



i r 00 ^ r 1 ,. f 2 \ , p d 6 



<t>f(x) = ^J^dydzJ dhdt 2 J 



|e — y||e — z||y — z| dx t 



x 2 



e x(y+z+\y-z\) _ i 



4 V 6 f e ip ' x f e iq ' x ^jC' 1 



7T^ 



t R R 2 (q-R) 2 ( P + R) 2 • (4 ' 16) 



2 T 6 J p p 2 7 q q 2 r R R 2 (q-R) 2 (p + R) 2 ' 



This corresponds to a contribution from graph (vii) of fig. Q] where a non-zero Matsubara 
mode circulates in one of the triangles; even if we shrink this "heavy-mode contribution" to 
a series of higher-order operators, there still remains an infrared contribution in which the 
remaining two propagators have Matsubara zero modes. These power-law contributions can 
be obtained by expanding the result of the sum-integral over R in powers of p and q and 
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Figure 3: The basic functions denned in eqs. (|4.1ip . (14. 12)) . together with their asymptotic limits as 
given m eqs. (j4T3| . (|4~T5l) . (|4~T7l) . The functions VY/m indicate the nature of the limit, and the 
argument shows the order up to which terms have been included. 



omitting contact terms, which immediately leads to 



^ R (x) 



360C(3) | 2520C(5) | Q 



x 



10 



X 



12 



X » 1 . 



(4.17) 



4.2. Asymptotic behaviours 



Inserting the asymptotic expansions in eqs. ([4^8]) . (|4.10[> . (|4.13l) . (|4.14j) . (|4.15p . (|4.17l) into 
eqs. (|4.5() . (I4.6h . we can write down the asymptotic expansions for the functions G^ and 
at small and large x. In the short-distance limit we get 
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(4.18) 



(4.19) 



Note that the dominant terms, i.e. the first rows, are temperature-independent. The ex- 
pressions obtained agree with the OPE results that can be extracted from ref. [12], but two 
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features are worth stressing. First of all, the leading thermal correction is of 0(g e ) and 
positive, whereas 0(<? 4 ) only gives a subleading thermal contribution at short distances; in 
equal-time correlators, in contrast, both orders contribute to the leading term 0{— 1/x 4 ), 
which is negative at short distances [61 [12] . This drastic difference illustrates that averaging 
over the Euclidean time coordinate has a significant effect if x < (5. Second, the OPE contri- 
butions proportional to (e — 3p)(T), which come with opposite signs in the two channels and 
distinguish them from each other [11} I12j . are of 0(g s /x s ) and thus beyond the accuracy of 
our present 0(g 6 ) analysis. 

At long distances, on the other hand, we obtain 
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(4.20) 



(4.21) 



The first rows indicate that at sufficiently large distances, x 2 ^> ir 2 /g 2 N c , the correction of 
0(g & ) overtakes the term of 0(g 4 ); this signals a breakdown of the perturbative series and 
necessitates a resummation, to which we now turn. 



5. Effective theory computation 



As can be observed from eqs. (|4.20|) . (|4.21|) . at distances of the order x 2 ~ \/g 2 T 2 N c the naive 
perturbative series breaks down and needs to be resummed in order to recover the correct 
long-distance behaviour. A systematic way to implement these resummations is provided by 
effective theories. In our case the relevant low-energy effective theory is known as EQCD, a 
three-dimensional theory for the Matsubara zero modes of Af and Aq [18\ [T9] . 

Keeping the dimensions of the fields as they are in the original theory, and continuing to 
denote by g 2 the renormalized 4d gauge coupling, the operators we are interested in can be 
represented in EQCD as 



ce9 



X E = c x g 2 



Z ¥ F^ + 2Z D (V t A ) a (V l A Q ) a + ... 
-4Z E e ljk {V,-M a Ff k + . . . 



(5.1) 
(5.2) 



where the numerical coefficients have been chosen so that the renormalization factors Z-p, 
Zj), i?E are all of the form 1 + 0{g 2 ). Dropping again contact terms, and setting e = since 
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the contributions are ultraviolet finite, the leading-order results can be written as 
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(47r) 2 x vdx 2 
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(5.4) 



where m\ is the Debye mass parameter appearing in the EQCD Lagrangian. To avoid double- 
counting, we write the observables as 



G e , x (x)= GlJx)-GtJx) 



+ 



GfJx) 



(5.5) 



where the quantity in the first square brackets is infrared safe. In the "naive" results, m E is 
treated as a perturbatively small quantity of O(g), and we can expand with respect to it in 
eqs. (15. 3p . (15. 4ft . giving 



( Gf(x) 
\±d A c 2 g± 

-lGdAclg 4 



7 2 
^F 



48vr 4 T 7 
x 6 

96tt 4 T 7 



48vr 4 T 7 



4vr 2 T 5 m 2 



4^ 2 T 5 



??? 



+ Ofm 4 



+ £>fm 



(5.6) 
(5.7) 



Substituting m E = g T N c /3, Zp = Zr> = Ze = 1 and comparing with eqs. (|4.20p and 
(|4.2ip . we see that these reproduce the dominant LO and NLO long-distance behaviours. 
Therefore, in both channels, the subtraction-addition step of eq. (|5.5p replaces the leading 
infrared contribution ~ —g e N c /x by an exponentially suppressed behaviour, proportional 
to ~ exp(— m E x) in the pseudoscalar (x) and to ~ exp(— 2m E x) in the scalar (9) channel, 
thereby postponing the breakdown of the perturbative series to larger distances. 

As far as the subleading infrared terms are concerned, we would expect those of 0(g & N c /x~ 5 ) 
to also be replaced by exponentially suppressed terms in the pseudoscalar channel (x), in 
which symmetries require the appearance of an odd number of Aq propagators in any in- 
termediate state [5]. (In other words, the effective theory for the static colour-magnetic 
modes, MQCD, does not contain operators with the right quantum numbers to contribute 
to G x .) In contrast, in the scalar channel, power-law terms of 0(g e N c /x 5 ) remain over, 
and imply that our computation becomes unreliable at the distance scale where this term 
overtakes the leading-order term, i.e. when C(l/x 6 ) ~ 0(g 2 N c /it 2 x 5 ). This happens at the 
non-perturbative colour-magnetic scale x ~ 7r / (g 2 N C T); the corresponding behaviour could 
in principle be determined through a non-perturbative analysis of MQCD. 
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As far as the remaining power-law terms go, those of 0(g®N c T 7 /x 6 ) could be accounted 
for by a choice of the 0(g 2 ) corrections to the renormalization factors Zp, Z-£>, Zp. Un- 
fortunately, our present results do not suffice to fix all of these unambigously, but further 
matching computations would be needed. Higher powers still, 0(g®T 7 /x 8 ) etc, correspond to 
higher-dimensional operators within EQCD or MQCD. Since the focus of the present paper 
is on intermediate distances, we have not carried out a next-to- leading order analysis of the 
correlators on the EQCD side. Rather, we will use leading-order resummation, as implied 
by eqs. (|5.3p - (|5.7p . to gauge the distance scale where infrared effects become important, and 
declare our results as unreliable above this scale. Explicitly, our resummation amounts to 
the additional contributions 

4o 2 jV c " 1 „4™6t^ / A2 \ 2 / „-2x 
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to be added to eqs. (|4.5 



| ff 4 m|r /d 2 1 

167r 2 x \dx 2 
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6. Numerical evaluation 

6.1. Choice of the renormalization scale 

Upon adding together the full theory and EQCD contributions to the correlators, eqs. fj4.5j) 
and (|5.8p for the 6 channel and eqs. (|4.6h and (15. 9h for the \ channel, we are now ready 
to numerically estimate the behaviour of the functions in various regimes. We do this for a 
variety of distance scales, ranging from x <C 1/ttT to a; > 1/nT, which in particular implies 
that we should choose an "optimal" renormalization scale that adapts to the physical quantity 
and distance regime of interest. At short distances, it is natural to determine the optimal 
scale from the UV limit of the correlators, demanding that the leading A-dependent term in 
the NLO part of the results vanishes. According to eqs. (|4.18p and (|4.19p this yields 

„103/165- 7E „13/165- 7E 

A£ = , A* = . (6.1) 

x x x 

At long distances, on the other hand, we fix the scale from the logarithmic terms in eqs. (|4.20p . 
(|4.2ip . which correspond to NLO corrections to the factors Zp, Zd, Zp of sec. [5] and are not 
otherwise cancelled by the resummation. This gives 

A^ = e" 7E 4vrT , A£ = e" 6/11 - 7E 4vrT . (6.2) 

(Remarkably, in both channels, the short-distance and long-distance scale choices coincide at 
x = e 10 3/i65-in2 ^ q.9334.) At intermediate distances, we make the simple choice of defining 
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Figure 4: The behaviour of the full 9 (left) and % (right) correlators, multiplied by x 6 , at T = 3T C . 
Thick lines indicate the scale choice A = A opt (eq. (|6.3[) ) , thin lines variations by a factor of 2 around 
this scale. Both correlators vary rapidly with x and are dominated by a vacuum-like behaviour up to 
distances x ~ a few. Our results become inaccurate at x > 10, cf. fig. [H 

an "optimal" scale through 



Subsequently we vary our A parameter by a factor of 2 around this scale. 

As far as the gauge coupling is concerned, we solve it from the 2-loop renormalization 
group equation, and define 



Where needed, the Debye mass parameter is evaluated from the leading-order expression 
m| = g 2 T 2 N c /3 with the same g 2 inserted as elsewhere. 

6.2. Plots 

We start by inspecting the behaviour of the full correlators over a wide distance range. This 
leads to the plots shown in figs. [H [5] where we display 
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Figure 5: Like fig. H but at T — 12T C 



multiplied by x in order to increase the resolution, on a log-log scale at the temperatures 
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i , 1, 2. A power-law behaviour 



can be clearly observed at short distances, whereas the dependence of the result on the 
renormalization scale as well as the magnitude of the NLO correction become visible at long 
distances, indicating that the perturbative series becomes less trustworthy there. 

In order to observe more clearly the fine structures of the correlators, we next consider the 
functions 
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Here the vacuum-like parts are defined through the first lines of eqs. (|4.18p . (|4.19p with a 
particular scale choice: 
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Although it would be tempting to interpret these as genuine vacuum results, this is not correct 
within our fixed-order computation: the genuine vacuum results become non-perturbative at 
a distance x ~ 0.1 fm and this is reflected by the fact that with the scale choices of eq. (16. ip 
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Figure 6: The short-distance behaviours of the "thermal" parts of the 9 (left) and \ (right) correlators, 
multiplied by x 3 , compared with the corresponding OPE limits (the second lines of eqs. (|4.18|) . (|4.19jl ). 
at T = 3T C . Thick lines indicate the scale choice A = A opt , thin lines variations by a factor of 2 
around this scale. (Numerical evaluation is challenging at x <C 1 due to significance loss: not only the 
"vacuum" parts l/x 7 but also "thermal" parts behaving as l/x 5 and l/x 4 cancel in AGg, x at x< 1, 
cf. eqs. I|4TTD|) . Km . KTB ). 



the perturbative corrections grow out of control at x > 0.1 A^g -1 [12]. In contrast, in the full 
finite-temperature results of eq. (|6.5p a scale choice like that in eq. (|6.3|) is well-motivated, and 
at sufficiently high temperatures the result is perturbative (apart from MQCD contributions) 
even at distances x ~ fm. The definitions of eqs. (|6.6p - (|6,9p therefore only serve purposes of 
illustration. In any case, in figs. [6H8] these functions are plotted on various distance scales. 

Considering first the UV limits of the correlators, in fig. [6] we display our results for AGq 
and AG X multiplied by x 3 , with x ranging from to 6, and compare the results with the UV- 
limits in eqs. (|4.18p . (|4.19p . The results indicate that OPE-expanded results are applicable 
only in the range x < 1, and that there is a visible difference between the two channels in the 
regime 1 < x < 5. 

Next, motivated by the results of ref. [6], we inspect intermediate distances, 3 < x < 12. 
The results are plotted in fig. [71 The results suggest that the two channels are rather close 
to each other in the distance range 5 < x < 10. 

Finally, fig. [8] shows the long-distance behaviours of AGo and AG X multiplied by x 5 (so 
that, within our approximation, they approach constants at large x). We compare the curves 
to the result of a calculation where no resummation of the soft modes was carried out. We 
observe that the —l/x 4 behaviour of the unresummed result becomes fully visible only when 
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Figure 7: The "thermal" parts of the 9 (left) and x (right) correlators at intermediate distances, at 
T = 3T C . Thick lines indicate the scale choice A = A opt , thin lines variations by a factor of 2 around 
this scale. Correlations are slightly more pronounced in the x channel. 



x > 10. Further higher order or non-perturbative corrections, which would eventually turn the 
behaviour in both channels into exponential decay [5], can be assumed to become significant 
in the same range where resummation starts to dominate the behaviour, i.e. x > 10 at this 
temperature. 



7. Conclusions 

The purpose of this paper has been to "interpolate" between previous studies of the spatial 
correlation functions related to the operators 6 oc Tr [F^F^] and x oc Tr [F^F^] within pure 
SU(3) gauge theory at a finite temperature T ~ a few x T c . We find that the "ultraviolet 
limit", studied in ref. [12] (following refs. [61 [20] and other recent works), in which the 
correlators are computed with Operator Product Expansion type methods, works reasonably 
only at distances 2irTx < 1 (cf. fig. [6]). At the same time the "infrared limit", studied in 
ref. [5] and many subsequent works, starts to govern the behaviour at distances 2irTx > 10 (cf. 
fig. [8]). In between, neither of these simplified approaches is viable, and our results represent 
the most precise weak-coupling expressions available to date. Whether they themselves are 
quantitatively accurate remains to be seen, once lattice measurements for the very same 
quantities, following the lines that have recently been published in ref. [6], become available. 
Among the more detailed motivations of our study was the question, to what extent do the 
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Figure 8: The long-distance behaviours of the 9 (left) and \ (right) correlators, multiplied by x 5 , 
at T = 3T C . The results are compared with the "naive" expressions which do not include the resum- 
mations of eqs. (|5.8p . (|5.9|) . Thick lines indicate the scale choice A = A op , thin lines variations by 
a factor of 2 around this scale. The "naive" results become negative at long distances, whereas the 
resummed results remain positive and behave as ~ 1 /x 5 within our approximation. (This would turn 
into an exponential decay by higher order EQCD corrections or by non-perturbative MQCD effects.) 



two correlators differ from each other? While we have indeed found some disparities between 
the two channels, in the direction that all features are more pronounced in the correlator 
corresponding to Tr [F^F^], the difference is not as striking as that seen in the lattice study 
of ref. Several possible reasons come to mind. First of all, the leading "conformality- 
breaking" OPE contribution proportional to (e — 3p)(T), showing a clear distinction between 
the two channels at short distances, corresponds actually to the 3-loop order in terms of 
Feynman graphs (in ref. [12] it could be "guessed" from the 2-loop computation carried out, 
by making use of renormalization scale invariance in connection with the general factorization 
philosophy inherent to OPE). This might suggest that 3-loop contributions are important 
insofar as the difference between the two channels is concerned. Second, it must be kept in 
mind that the lattice measurements in ref. [6] concerned strictly local correlators, whereas 
ours have been averaged over the Euclidean time coordinate. This has an important effect 
at least at distances x</3, i.e. 2ttTx < a few. It will be interesting to repeat the comparison 
if lattice results for time-averaged correlators become available. Third, it should be stressed 
than in ref. [6] the genuine vacuum parts were subtracted from the thermal correlators. 
Unfortunately, this makes the difference sensitive to confinement physics; the features seen 
may not be a reflection of perturbative thermal physics alone. It would be useful to carry 
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out future comparisons without such a subtraction, i.e. directly with figs. 01 

Let us remark that the "thermal" parts of the correlators, as well as the difference be- 
tween the two "thermal" parts, show a remarkably rich x-dependence, suggesting that a high 
resolution is needed for disentangling all relevant features. As can be observed in figs. [6HE1 
the "thermal" parts change sign twice, before finally showing a positive correlation at long 
distances. At the same time, the difference between the two "thermal" parts, proportional to 
the thermal part of the function </>lo (cf. eq. (I4.6P and fig. [2]), shows the same two sign changes 
at short and intermediate distances, but would have a third one in store for long distances: 
at the "realistic" temperatures considered, the screening mass corresponding to Tr \F^ il jF^ il j\ 
is markedly larger than that corresponding to Tr [F^F^] (cf. refs. [HI E]), implying that 
at very large distances the correlator G x would be smaller in magnitude than Gq. 

We end by recalling the relation of our work to heavy quarkonium physics and simultane- 
ously the physical reason for considering time-averaged correlators as we have done in the 
present study. Whereas it was traditionally thought that heavy quarkonium melting is pri- 
marily a consequence of colour-electric Debye screening at long distances, x ~ it has 
been suggested in recent years that at least in the bottomonium case, where weak-coupling 
techniques are expected to be semi-quantitatively applicable, the "melting" of quarkonium 
resonances takes place already when the Bohr radius of the bound state, r, is parametrically 
in the range 1/ttT C r < |21|-|24j. In this situation the dominant mechanism for 
melting is not Debye screening but Landau damping or, more physically, the decoherence of 
the quantum mechanical bound state, mathematically represented by an imaginary part in 
the real-time static potential felt by a heavy quark-antiquark pair propagating in Minkowski 
time [25] (see also refs. [26]-[29]). Moreover, it can be argued that it is interesting to study 
quarkonium states even before they melt, whence the regime x<1/ttT becomes relevant 
as well [30]. With such motivations, the "singlet free energy" as determined from lattice 
simulations [31] and traditionally used in potential models has recently been compared with 
perturbation theory at intermediate distances x ~ 1/vrT, and a surprising agreement between 
the two was found |32] Q The singlet free energy consists precisely of a correlator of two time- 
integrated objects. It is for this reason that we suspect time-averaged correlators to be of 
more direct physical relevance than strictly local ones. By learning about the reliability of 
the weak-coupling expansion for observables of this type in the intermediate distance range, 
it will hopefully become possible in the end to relate perturbative determinations of heavy 
quarkonium spectral functions (ref. [33] and references therein), which use the real-time static 
potential as input, to on-going lattice studies such as those in refs. [351 [36], thereby learning 
about the systematic errors on both sides. 

A similar if inequivalent observable, the correlator of two "traced" Polyakov loops, sometimes referred 
to as (exponent of minus) the "colour-averaged potential", was computed perturbatively in ref. [33], but no 
numerical evaluation or comparison with lattice data was presented. 
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